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1.  INTRODUCTION 


1.1  FOCUS  OF  THE  PILOT  STUDY 

Short  range  forecasts  for  periods  of  24  hours  or  less 
usually  are  based  on  both  the  objective  guidance  provided  by 
large  scale  models  and  the  individual  forecaster's  experience 
gained  over  time  at  a  particular  location.  The  forecaster  who 
lacks  this  personal  experience,  however,  often  must  rely  on 
"textbook"  interpretation  of  a  model's  objective  guidance,  and 
thus  possibly  be  handicapped  in  forecasting  primarily  mesoscale 
si tuati ons . 

This  study  examines  the  skill,  usefulness  and  applicability 
of  two  techniques  that  have  been  applied  to  short  range  forecasting 
as  substitutes  for  local  experience.  They  are  known  as  regression 
estimates  of  event  probability  (REEP),  and  logistic  regression  (LR) 
probabi 1 i ties . 

These  techniques  are  based  on  conditional  climatologies  (see 
Miller,  1979).  They  use  surface  observations  coupled  with 
forecast  equations  based  on  regression  techniques,  to  produce 
probability  forecasts. 

The  skills  of  the  techniques  in  forecasting  summer  stratus 
ceilings  at  four  naval  air  stations  in  California  --  Moffett  NAS, 
North  Island  NAS,  Miramar  NAS,  and  San  Nicolas  NAS  --  were 
evaluated  in  this  study.  Both  REEP  and  LR  were  used  at  Moffett; 
only  LR  was  used  at  the  other  three  stations.  (REEP  equations 
also  were  used  operationally  at  Moffett  NAS  during  the  summer  of 
1982;  see  Appendix  A  for  results  of  this  evaluation.)  Skill 
scores  of  each  technique  were  compared  to  those  of  climatology,  a 
readily  quantifiable  technique,  to  assess  performance. 
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1.2  THE  FORECAST  PROBLEM 

Summer  stratus  along  the  California  coast  presents  a  daily 
forecast  problem.  Depending  on  the  location,  stratus  with 
ceilings  as  low  as  500  ft  can  develop  during  the  late  afternoon 
or  evening,  dissipating  sometime  the  following  morning.  A  series 
of  days  with  stratus  occurring  typically  is  followed  by  one  to 
several  clear  days;  California  summer  is  typically  either  coastal 
stratus  or  clear. 

Since  the  mid-latitude  Pacific  anticyclone  is  a  persistent 
feature  in  the  summer,  and  the  large  scale  synoptic  pattern  does 
not  change  much,  the  forecast  problem  therefore  is  associated 
with  subtle  changes  in  the  mesoscale  features.  Predictors  such 
as  surface  winds,  dew  point  depression,  etc.,  are  strong 
candl dates . 

Topographical  features  also  play  a  major  role  since  the 
stratus  deck  is  confined  to  the  marine  layer.  Cool  sea  surface 
temperatures  along  the  California  coast  in  summer  cause  a  low 
level  temperature  inversion  and  an  associated  moist  marine  layer. 
The  predominant  on-shore  pressure  gradient  forces  this  marine 
layer  on  shore,  with  diurnal  modification  from  the  sea  breeze 
gradient.  Coastal  terrain  features  modify,  or  dam  up  this  inland 
intrusion. 

Moffett  NAS,  for  example,  is  located  on  the  southern  edge  of 
the  San  Francisco  Bay,  with  the  Santa  Cruz  range  to  the  west 
(Figure  1).  These  hills  cause  the  stratus  to  approach  Moffett 
from  the  north,  along  the  bay.  The  cool  and  moist  marine  layer 
must  slide  through  the  Golden  Gate  and  push  down  the  bay,  so  the 
stratus  deck  that  occurs  at  Moffett  is  really  a  combination 
advection-formation  event. 
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Figure  1.  Locations  of  air  stations  and  airports,  with 
mountain  elevations  to  the  west  and  east,  in  the  San 
Francisco  Bay  Region  of  Northern  California. 
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At  the  three  air  stations  in  Southern  California,  the 
stratus  deck  has  a  clear  path  inland;  field  elevation  and 
distance  inland  are  the  primary  differences  (see  Figure  2).  The 
lower  latitude  and  much  warmer  sea  surface  temperatures  along  the 
southern  coast  also  have  an  effect,  causing  thicker  stratus  decks 
with  higher  bases. 


2.  DESCRIPTIONS  OF  THE  REEP  AND  LR  TECHNIQUES 

2.1  REGRESSION  ESTIMATES  OF  EVENT  PROBABILITIES  (REEP) 

The  REEP  technique  (described  in  detail  in  Miller,  1964) 
uses  surface  observations  from  the  forecast  station  and  other 
surrounding  stations.  Each  predictor  variable  is  divided  into 
subranges;  these  are  then  assigned  values  of  0  or  1,  depending  on 
the  observed  values. 

(Lists  of  candidate  predictors  for  both  REEP  and  LR 
(logistic  regression)  techniques  used  in  this  pilot  study  are 
given  for  Moffett,  North  Island,  Miramar  and  San  Nicolas  Naval 
Air  Stations  in  Appendixes  A-E  respectively.) 

As  an  example,  the  predictor  variable  "surface  dew  point 
depression"  is  divided  into  five  subranges  for  Moffett  NAS  data: 

Dummy  Variable  No.  Range  (°F) 

1  -8.0 

2  8.1-10.0 

3  10.1-12.0 

4  12.1-14.0 

5  >14.1 

The  dew  point  depression  variable  can  now  be  represented  by 
a  vector, 

D  =  (1,0, 0,0,0), 

which  indicates  the  value  is  in  the  0-8. 0°F  range.  All  the  other 
predictor  variables  are  categorized  similarly,  and  the  resulting 
single  surface  observation  for  any  given  day  and  time  becomes 
simply  a  string  of  0's  and  l's. 


Using  linear  regression  methods,  forecast  indices  are  found 
for  each  forecast  category  from  the  regression  equations,  one  for 
each  forecast  category: 


=  a 


Oi 


+  a 


li  *1 


+  \i  V  1  =  ’•  2 


The  F.j's  are  the  forecast  indices  for  each  ceiling  category: 
Fi  for  ceilings  less  than  or  equal  to  1000  ft  and  F2  'for  ceilings 
greater  than  1000  ft,  for  example.  The  a^  represent  the  regres¬ 
sion  coefficients  and  the  xi  represent  the  0,1  dummy  variables. 
Although  more  than  two  ceiling  categories  can  be  selected, 
greater  skill  results  when  only  two  are  used. 

Fi  resembles  a  probability  in  two  ways:  F i  +  F 2  =  1;  and  a 
larger  Fi  means  that  state  i  is  more  likely.  However  individu  ' 

Fi  are  not  bounded  by  0  and  1.  Therefore  Fi  is  not  rigorousl, 
probability  and  will  be  referred  to  as  a  forecast  index. 

An  additional  characteristic  of  REEP  is  that  since  the  x 

equal  0  or  1,  the  magnitude  of  the  a  -j  j  directly  show  the  pred. 
tors  relative  influence  on  the  index  value. 


2.2  LOGISTIC  REGRESSION  (LR) 

Consideration  of  a  forecast  parameter  which  has  the 
properties  of  a  probability  leads  to  the  following  (logistic) 
transf ormati on ; 


where 


and 


N 

A .  =  l  a.-  x.;  N  =  No.  of  predictors 
1  j=l  1J  J 


0  . 


In  (ttB  -  X, 


is  called  the  log-odds  ratio. 
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1,  therefore  6  satisfies  the 


Now  CK9 ^£l  and  0^+02  = 
probability  properties. 

The  coefficients  a^j  are  found  from  applying  regression 
techniques  to  the  dependent  data  sample  (see  Dixon,  1981).  The 
advantage  of  LR  is  that  continuous  variables  can  be  used  without 
first  categori zi ng;  also,  categorical  variables  can  be  used 
directly  without  transformation.  However,  only  binomial  predic- 
tands  are  al 1  owed. 

Predictor  selection  is  based  on  a  statistical  stepwise 
procedure  called  asymptotic  covariance  estimate.  An  explanation 
of  this  procedure  is  beyond  the  scope  of  this  report,  and  the 
reader  is  referred  to  Dixon,  1981,  for  information. 

3.  RESULTS 

3.1  COMPARISONS  OF  REEP  AND  LR  AT  MOFFETT  NAS 

The  REEP  and  LR  equations  are  given  in  Appendices  A  and  B, 
respectively,  for  ceilings  1000  ft  and  2000  ft. 

Five  forecast  valid  times  (VT)  were  selected  for  this  pilot 
study;  a  common  forecast  generation  time  of  1800  LST  allowed  use 
of  the  0000  GMT  Oakland  radiosonde  data,  and  also  was  a  time 
prior  to  stratus  formation.  The  equations  for  REEP  were 
generated  from  ten  years  of  data;  the  equations  for  LR  were 
generated  from  seven  years  of  data. 

Tables  1  and  2  show  percent  correct  and  threat  scores  for 
REEP,  LR,  and  climatology  for  the  dependent  and  independent  data 
samples,  respectively. 

Threat  score  is  defined  as 

TS  =  A/(A+B+C)  where 

A  is  the  number  of  correct  stratus  forecasts,  B  is  the 
number  of  incorrect  non-stratus  forecasts,  and  C  is  the  number  of 
incorrect  stratus  forecasts. 
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Table 

1 . 

Dependent  data 

compari sons . 

Cei  1 

ing  <2000 

Ft 

Percent 

Correct 

Threat  Score 

VT 

REEP 

LR 

CLIMO 

REEP 

LR 

CLIMO 

2200 

88 

83 

88 

.34 

.37 

.06 

2400 

77 

83 

72 

.24 

.55* 

.16 

0200 

77 

82 

67 

.57 

.59 

.20 

0400 

75 

78 

50 

.61 

.57 

.28 

0600 

78 

77 

56 

.64 

.58 

.28 

Ceil 

ing  _<1000 

Ft 

Percent 

Correct 

Threat  Score 

VT 

REEP 

LR 

CLIMO 

REEP 

LR 

CLIMO 

2200 

92 

82 

92 

.25 

.23 

.04 

2400 

85 

85 

85 

.32 

.46* 

.08 

0200 

84 

84 

77 

.51 

.55 

.13 

0400 

81 

80 

73 

.48 

.51 

.16 

0600 

78 

77 

72 

.47 

.46 

.16 

8 


Table  2.  Independent  data  comparisons  (200  observations) 


Cei  1  i  ng  £2000  Ft 


Perc 

ent 

Cor rec  t 

Threat  S 

core 

VT 

REEP 

LR 

CLIMO 

REEP 

LR 

CLI 

MO 

2200 

88 

82 

8  7; 

.34 

.37 

.07 

2400 

78 

82 

75 

.27 

.57* 

.14 

0200 

76 

80 

63 

.54 

.60 

.22 

0400 

75 

68 

57 

.62 

.42 

.27 

0600 

74 

76 

52 

.57 

.62 

.32 

Cei 

1  i ng  £1000  F 

t 

Perc 

ent 

Correct 

Threat  Score 

VT 

REEP 

LR 

CLIMO 

REEP 

LR 

CLI 

MO 

2200 

91 

90 

92 

.24 

.09* 

.04 

2400 

86 

84 

86 

.24 

.47* 

.07 

0200 

72 

78 

76 

.38 

.50 

.14 

0400 

75 

62 

70 

.38 

.43 

.18 

0600 

70 

70 

69 

.38 

.47 

.18 

r* 
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Given  this  relationship,  0£TS£1  with  TS  =  1  represents  a 
perfect  score  while  TS  =  0  indicates  a  no-skill  score. 

For  ceilings  £2000  the  percent  correct  comparisons  show  both 
REEP  and  LR  are  better  than  climatology  for  the  last  four  valid 
times  (VT).  LR  also  seems  to  be  slightly  better  than  REEP.  For 
ceilings  £1000  ft,  however,  no  clear  advantage  can  be  seen  for 
either  method. 

The  threat  score  comparisons  show  both  REEP  and  LR  have 
definite  skill  over  climatology*,  with  the  exception  of  the  2200 
VT  LR  in  Table  2b.  A  contingency  table  for  2200  VT  LR  can  be 
shown  as 


Forecast 

Ceiling  Ceiling 


Observation 

£1000  Ft 

>1000  Ft 

Total 

Ceiling 

£1000 

Ft 

2 

14 

16 

Cei 1 i ng 

>1000 

Ft 

6 

178 

184 

Total 

8 

192 

200 

With  the  exception  of  the  LR  threat  score  for  the  2200  VT  in 
Table  2b,  the  only  LR  scores  significantly  better  than  REEP  (5% 
chi  square  test)  are  for  the  2400  VT  forecast.  This  exception  is 
easily  explained:  the  above  tabulation  shows  a  small  number  of 
observed  ceilings  £1000  ft;  therefore,  a  small  fluctuation  in 
the  number  of  correct  low  ceiling  forecasts  can  account  for  the 
low  threat  score  for  the  2200  VT  LR  forecast  in  Table  2b. 


*Climatology  threat  score  =  P / ( 2 - P )  where  P  is  the  climatological 
probabi 1 i ty . 


The  explanation  for  the  universally  better  2400  VT  LR  threat 
scores  is  not  so  clear.  With  the  exception  of  Table  2b,  the  REEP 
2400  VT  scores  seem  inconsistent  with  the  other  REEP  threat 
scores.  As  the  valid  time  increases,  the  frequency  of  low 
ceilings  increases.  The  threat  scores  characteristically  become 
larger  with  increasing  relative  frequency.  This  trend  in  REEP 
threat  score  is  not  followed  at  the  2400  VT,  even  thouyh  the 
relative  frequencies  are  nearly  double  those  of  2200  VT. 

It  must  be  assumed  that  the  data  used  for  REEP  equations  had 
some  anomalous  characteristic  at  the  2400  observation;  the  2400 
VT  equations  therefore  will  not  be  considered  in  this  comparison. 

The  following  conclusions  can  be  drawn  for  the  REEP/LR 
comparison  at  Moffett  NAS: 

(1)  Percent  correct  scores  snow  that  both  REEP  and  LR  are 
better  than  climatology,  except  for  the  2200  VT  forecasts. 

(2)  Threat  scores  are  clearly  better  for  both  REEP  and  LR 
as  compared  to  climatology. 

(3)  After  the  2400  VT  forecast  is  eliminated,  REEP  and  LR 
show  about  the  same  skill  scores. 

3.2  REEP  OPERATIONAL  EVALUATION  AT  MOFFETT  NAS 

REEP  operational  evaluation  results  are  given  in  Appendix  A. 
REEP  was  used  since  Moffett  did  not  have  the  LR  program  available 
at  the  time  of  the  study.  The  subjective  conclusions  of  the 
evaluation  were  (1)  Reep  was  easy  to  use;  (2)  REEP  has  a  skill 
about  equal  to  that  of  an  experienced  forecaster;  and  (3)  REEP 
does  not  improve  a  good  forecaster's  skill.  (Close  data  examina¬ 
tion  reveals  inconsistencies  which  cannot  be  completely  resolved 
with  the  existing  data;  see  Appendix  1.) 


11 


3.3  LR  RESULTS  FOR  NORTH  ISLAND,  MIRAMAR  AND  SAN  NICOLAS  NAS 


Appendices  C,  D,  E  give  the  LR  equations  for  North  Island, 
Miramar,  and  San  Nicolas  NAS. 

Tables  3  through  5  give  the  dependent,  independent  and 
climatology  skill  scores  for  the  three  NAS. 

Table  3a.  North  Island  -  ceiling  £3500  Ft. 

Percent  Correct  Threat  Score 


LST 

DEP 

IND 

CLIMO 

DEP 

IND 

CLIMO 

Relative 

F  requency 

1900 

67 

72 

62 

.49 

.52 

.23 

.38 

0100 

75 

79 

67 

.73 

.76 

.50 

.67 

0700 

84 

74 

76 

.83 

.74 

.61 

.76 

Table  3a  shows  that  (1)  threat  scores  are  better  than 
climatology;  (2)  threat  scores  increase  with  lead  time;  (3) 
percent  correct  is  not  significantly  better  than  climatology;  and 
(4)  relative  frequency  increases  with  lead  time. 

Table  3b  reveals  the  same  information  as  Table  3a,  with  the 
exception  of  point  3).  In  fact,  percent  correct  overall  is  not 
as  high  as  climatology,  which  clearly  diminishes  the  forecast's 
worth  despite  of  the  better  threat  scores. 


Table  3b. 

North 

Island  - 

ceiling 

<1000 

Ft. 

Percent 

Correct 

Threat 

Score 

LST 

DEP 

IND 

CLIMO 

DEP 

IND 

CLIMO 

Rel ati ve 

F  requency 

1900 

90 

80 

93 

.29 

.10 

.04 

.07 

0100 

73 

57 

76 

.38 

.34 

.14 

.24 

0700 

69 

64 

70 

.43 

.36 

.18 

.30 

Tables  4a  and  4b  show  that  (1)  threat  scores  are  better  than 
climatology;  (2)  percent  correct  about  the  same  as  climatology; 
(3)  relative  frequencies  for  the  0100  LT  and  0700  LT  forecasts 
are  about  the  same;  and  (4)  threat  scores  for  0100  LT  and  0700  LT 
forecast  are  about  the  same. 

Table  4a.  Miramar  -  ceiling  £3500  ft. 


Percent 

Correct 

Threat 

Score 

LST 

DEP 

IND 

CLIMO 

DEP 

IND  CLIMO 

Rel ati ve 

F  requency 

1900 

76 

82 

72 

.46 

.58 

.16 

.28 

0100 

76 

76 

64 

.71 

.72 

.47 

.64 

0700 

79 

75 

69 

.77 

.72 

.53 

.69 

Table  4b . 

Mi ramar 

-  ceiling  £1000 

ft. 

Percent 

Correct 

Threat 

Score 

LST 

DEP 

IND 

CLIMO 

DEP 

IND  CLIMO 

Rel ati ve 

F  requency 

1900 

86 

90 

92 

.17 

.26 

.04 

.08 

0100 

66 

65 

61 

.51 

.52 

.24 

.39 

0700 

61 

56 

62 

.50 

.47 

.23 

.38 

I  n 

Tables 

4a  and 

4b  the  scores  for 

the  LR 

again 

shows 

considerable  skill  over  climatology.  However,  the  similar 
percent  correct  for  the  two  LR  samples  and  climatology  points  out 
a  weak  overall  forecast  skill. 

Tables  5a  and  5b  show  that  (1)  threat  scores  are  consistent¬ 
ly  higher  than  climatology;  (2)  percent  correct  is  better  than 
climatology,  but  lower  than  for  the  other  three  NAS;  and  (3) 
relative  frequencies  are  higher  for  San  Nicolas  than  for  North 
Island  or  Miramar. 
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Table  5a. 

San 

Nicolas  - 

ceiling 

<3500 

ft. 

Percent 

Correct 

Threat 

Score 

LST 

DEP 

IND 

CLIM0 

DEP 

IND 

CLIMO 

Relative 

F  requency 

1900 

73 

71 

61 

.69 

CO 

LO 

• 

.43 

.61 

0100 

76 

73 

68 

.73 

.70 

.51 

.68 

0700 

81 

80 

76 

o 

CO 

• 

.80 

.61 

.76 

Tabl e  5b. 

San 

Nicolas  - 

ceiling 

<1000 

ft. 

Percent 

Correct 

Threat 

Score 

LST 

DEP 

IND 

CLIMO 

DEP 

IND 

CLIMO 

Relative 

F  requency 

1900 

69 

69 

58 

.52 

.59 

.27 

.42 

0100 

68 

68 

53 

.56 

.60 

.31 

.47 

0700 

67 

64 

52 

.55 

.55 

.32 

.48 

Because  of 

point 

(3)  relative  to 

Tables  1 

5a  and 

5b,  very 

little  variation  occurs  between  the  dependent  and  independent 
samples  scores.  This  fact  suggests  that  the  equation's  skill 
merely  reflect  the  underlying  climatology,  and  that  differences 
between  dependent  and  independent  samples  imply  climatological 


differences  in  the  samples  themselves. 

4.  CONCLUSIONS  AND  RECOMMENDATIONS 

Short-range  stratus  ceiling  forecast  equations  were 
generated  at  four  naval  air  stations  on  the  California  coast 
using  logistic  regression  (LR)  probabilities.  A  second  statis¬ 
tical  technique  --  regression  estimates  of  event  probabilities 
(RF.EP)  --  was  used  to  develop  Moffett  NAS  forecast  equations. 
These  equations  were  evaluated  operationally  at  Moffett;  when  the 
logistic  program  became  available  later,  they  were  compared 
to  LR. 


4.1  CONCLUSIONS 

The  REEP  operational  evaluation  at  Moffett  indicated  that 
although  the  technique  has  definite  skill,  the  experienced 
forecaster  could  do  as  well  not  using  it  as  an  inexperienced 
forecaster  could  do  using  it.  The  concensus  at  Moffett  was  that 
REEP  would  be  helpful  to  new  duty  forecasters,  but  that  they 
would  place  less  reliance  on  REEP  after  they  gained  more 
experience.  The  primary  shortfall  of  both  REEP  and  experienced 
forecasters  lay  in  forecasting  the  transition  from  stratus  to 
non-stratus  sequence,  and  vice  versa. 

The  REEP-LR  comparison  at  Moffett  showed  essentially  no 
difference  in  the  two  skill  scores.  Both  were  significantly 
better  than  climatology  in  the  th  aat  score  comparisons,  but  only 
marginally  better  in  most  percent-correct  comparisons. 

The  LR  equations  for  the  three  Southern  California 
stations  --  given  in  Appendices  C,  D  and  E  --  yield  unimpressive 
results.  The  threat  scores  were  universally  better  than  climato¬ 
logy,  but  the  percent-correct  scores  were  marginal,  with  some 
better  and  some  worse  than  climatology. 

Overall,  the  study's  results  indicate  that  the  statistical 
forecast  techniques  are  only  marginally  useful  as  forecasting 
tools  in  the  particular  applications  investigated  here. 

Certainly  climatology  can  not  be  viewed  as  wholly  valid 
competition,  insofar  as  the  actual  competition  for  any  forecast 
procedure  would  be  the  operational  forecaster's  skill  scores 
(none  available  for  this  pilot  study). 

The  stations  do  not  have  a  universal  forecast  verification 
procedure;  each  office  sets  its  own,  and  apparently  can  revise  it 
at  will.  It  becomes  very  difficult,  therefore,  to  test  a 
forecast  procedure  for  operational  usefulness  when  no  long-term 
and  rigorous  skill  scores  exist  for  comparison. 
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At  locations  such  as  foreign  coastal  regions  where  there  is 
archived  data  but  no  U.S.  forecast  experience,  techniques  like 
REEP  or  LR  obviously  would  be  useful.  Two  criteria  must  be  met, 
however:  (1)  the  location  must  have  about  10  years  of  data 
available,  and  (2)  the  selected  forecast  parameters  must  have  a 
significant  dependency  on  local  surface  conditions.  If  the 
desired  forecast  parameter  were  rainfall  intensity,  then  clearly 
these  procedures  would  be  useless.  Parameters  such  as  surface 
winds  and  visibility,  however,  would  conceivably  be  candidates 
for  such  autonomous  statistical  procedures. 


4.2  RECOMMENDATIONS 

Based  on  the  information  and  experience  developed  during  the 
pilot  study,  the  following  are  recommended: 

(1)  Establish  universal  forecast  verification  procedures  at 
each  NAS  forecast  office.  These  procedures  should  not  be  based 
only  on  simple  percent  correct  (since  such  scores  can  be  very 
misleading  by  themselves),  but  upon  threat  scores,  Heidke  skill 
scores,  and  meaningful  forecast  parameters  such  as  ceiling 

hei ght . 

(2)  Compile  and  archive  verification  statistics  so  that 
prospective  forecast  techniques  can  be  compared  prior  to  opera¬ 
tional  evaluation.  The  operational  evaluation  should  be 
performed  mainly  to  iron  out  any  implementation  problems  and 
special  situation  difficulties. 

(3)  After  a  period  of  time,  conduct  a  survey  to  determine 
areas  of  potential  improvement  in  forecasts. 
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APPENDIX  A 


EVALUATION  OF  REEP  EQUATIONS  AT  MOFFETT  NAS 
A . 1  RESULTS 

Ceiling  forecasts  for  2200,  2400,  0200,  0400  and  0600  local 
time  (LT)  were  developed  for  Moffett  using  regression  estimates 
of  event  probabilities  (REEP).  The  complete  set  of  equations  is 
given  in  Section  A. 2.  Forecasts  for  ceilings  above  or  below  1000 
ft  are  generated  using  1800  LT  observations  and  Oakland  0000  GMT 
(1600  LT)  radiosonde  data.  The  ceiling  forecasts  are  limited  to 
the  June-September  period  of  summer  coastal  stratus  occurrence  at 
this  location. 

The  evaluation  began  in  June  1982  and  used  the  forecast 
packet  reproduced  as  Section  A. 3.  Beginning  in  July  with  the 
arrival  of  a  HP-41CV  calculator  and  printer,  the  forecast  proce¬ 
dure  was  automated  and  run  on  the  calculator.  An  output  example 
appears  in  Figure  A-l. 

Each  day  the  duty  forecaster  recorded  these  forecasts  and 
verifications  for  the  previous  day.  Also  noted  was  whether  REEP 
was  used  alone  or  with  additional  guidance,  or  the  forecast  was 
based  solely  on  other  guidance.  Tables  A-l  and  A-2  provide  these 
record  s . 


Table  A-l.  Moffett  stratus  forecast. 


Forecast 

REEP 

REEP  and 

Made 

Alone 

Other 

Other 

A1 1 

99% 

83% 

57% 

Stratus 

★ 

46% 

44% 

♦None  observed 
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NnF0GT  FGRt'.  «  r  • 

FOP  2800"  ::  0.04 
VERIFY  PERCENT  1! 

0600GHT  F0?U'5I'T 
for  1060FT  ;e  e.es 

VERIFY  PERCENT  IE  9t. 
0300GHT  FOF'EEhS7 

for  imr.  ie  -0.es 

VERIFY  PERCENT  IE  92. 

0800GK’  FORECAST 
for  1000ft  ;c  e.e- 
VERIFY  PEFCENT  13  97. 

1800GHT  FORECAST 
FOR  2000F7  IS  t.2\ 
VERIFY  PERCENT  IS  89. 

1800GWT  FORECAST 
r0R  1000FT  I:  e.<: 
VERIFY  PERCENT  !':  9S. 

1200GPT  FORECAST 
C0R  2000"  IE  0. SE¬ 
VER  I FY  PEF.EEN-  IE 

:Z00GIC  FORETf :•* 

-OF  I000"T  IC  -8.04 
VERIFY  PERCENT  l-  93. 

1 4005PT  FOPT^I* 

FO?  2000F"  I:  0.29 
VERIFY  PERCE*’  IS  74. 

!400GN*  FOPZrc-* 

FOR  I000FT  tc  0,72 
VERIFY  PERCENT  IS  94. 


Figure  A-l.  HP-41CV  output  for  13  August  1  982  giving  i£2000  ft 
and  -£1000  ft  forecasts  for  all  five  forecast  times. 
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Table  A-2.  Percent  correct  forecasts  (forecaster 
"agreed"  or  "disagreed"  with  REEP). 

Forecast 

Made  Agreed  Disagreed 

Non-stratus  99.5%  * 

Stratus  83%  24% 

*None  observed 

Table  A-l  gives  the  verification  summary  for  the  two  types 
of  forecasts.  For  "all"  forecasts,  REEP  apparently  was  superior. 
For  the  non-stratus,  no  clear  advantage  resulted.  The  low  score 
for  the  "stratus"  forecasts  compared  to  the  "all"  forecasts 
results  from  the  difficulty  of  forecasting  the  stratus  event. 

Long  periods  of  clear  weather  were  followed  by  short  periods  of 
stratus,  and  the  initiation  of  such  stratus  periods  is  difficult 
to  forecast. 

Due  to  the  long-running  clear  sequence,  the  persistence 
factor  tends  to  influence  the  forecaster.  During  obvious  clear 
spells,  the  forecaster  is  more  likely  to  go  with  REEP  alone. 

There  also  seemed  to  be  some  forecasters  during  the  study  who 
never  trusted  REEP  and  always  indicated  that  they  used  additional 
guidance.  This  could  easily  contribute  to  the  8 3 % - 5 7 %  advantage 
over  the  "other"  category. 

The  results  in  Table  A-2  are  a  little  harder  to  explain  in 
the  light  of  the  unenthusi asti c  response  related  verbally  by 
forecasters.  The  somewhat  small  sample  size  may  have  had  some 
effect,  but  a  complete  explanation  is  unknown. 
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A. l  MOFFETT  REEP  EQUATIONS 

The  equations  are  functions  of  binary  predictor  variables 
Xi,  whose  values  are  “1"  if  the  variable  is  in  the  given  range 
and  "0"  otherwise.  The  are  listed  at  the  end  of  this 
secti on . 

Forecast  Equations  -  Veiling  <1000' 

Val  id  Time 

2200  LT 

F  =  . 469+. 183X.-.254X15_. 141X4 
2400  LT 

F  =  . 245+. 325X. - .245X4+. 107X19 

+  .096X20+ -075X23- • 146X15- • 104X9 

0200  LT 

F  =  . 379+.  352X.-.272X6+.  173X20--  141X24-. 099X22 
0400  LT 

F  =  . 325+. 414X.-.240X7+. 230X20 

-  .257X24-  .095X22-  -153X18  -  •  121X26 

-  .248X25 

0600  LT 

F  =  . 476-. 257X8+.  273X.-. 071X22 
+  .212X20- • 167X24 
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Forecast  Equations  -  Ceiling  £2000' 

2200  LT 

F  =  .  542+.  301X.-.  348X16_. 155x10 
2400  LT 

F  =  .165+.508X.-.338X16  +  . 220X2 

-.129XH+. 110X17+. 131X19 

0200  LT 

F  =  • 224+. 694X.  +  .212X2o-.200Xi2+. 185X2 

0400  LT 

F  =  . 31 1 +. 544 X .+. 286 X2q+ . 248X2 

- .130X13- . 121X3+.089X21+. 116X19 

0600  LT 

F  =  .  3 1 6+ . 685 X .  +  .  1 86 X2g - . 1 62 X 1 4 
-. 100X13-. 124X22+. 132X2 
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List  of  Predictors  (1800  LT  observation  unless  specified) 
At  Moffett 


X1 

& 

Dew  point 

depression  (0 

-8  °F ) 

X2 

= 

Dew  point 

depression  (8 

-10°F) 

X3 

Dew  point 

depression  (H4°F) 

X4 

= 

Previous 

Day' s 

Ceiling 

(>1000' ) 

at 

2200 

LT 

X5 

= 

Previ ous 

Day's 

Cei 1 i ng 

(>1000'  ) 

at 

2400 

LT 

X6 

= 

Previ ous 

Day' s 

Cei 1 i ng 

(>1000* ) 

at 

0200 

LT 

X7 

= 

Previ ous 

Day's 

Ceiling 

(>.1000'  ) 

at 

0400 

LT 

X8 

* 

Previ ous 

Day' s 

Ceiling 

(>1000'  ) 

at 

0600 

LT 

X9 

= 

Wind  Speed  (0-6 

i  kt) 

X 

o 

= 

Prev i ous 

Day' s 

Cei 1 i ng 

(>2000°) 

at 

2200 

LT 

X11 

= 

Previ ous 

Day' s 

Ceiling 

(>.2000°) 

at 

2400 

LT 

X 1 2 

3 

Prev i ous 

Day '  s 

Ceiling 

(>2000°) 

at 

0200 

LT 

X13 

8 

Previ ous 

Day' s 

Ceiling 

(>2000°) 

at 

0400 

LT 

X 1 4 

3 

Previ ous 

Day' s 

Ceiling 

(>2000°) 

at 

0600 

LT 

At  Alameda 


X15  3 

Cei  1  i  ng  ( M0QQ  1  ) 

X 1 6  * 

Ceiling  (>2000') 

X 1  7  3 

Dew  Point  Depression 

(0-8°F ) 

X18  s 

Dew  Point  Depression 

(>14°F) 

X19  3 

Wind  Speed  (12-15  kt) 

At  San  Francisco 

*20  3  Dew  Point  Depression  (0-8°F) 
X 2 1  *  Wind  Direction  (  282-326) 

*22  3  Wind  Direction  (237-281) 
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Oakland  Radiosonde  (0000  GMT) 

*23  *  Temperature,  950  mb  (20-25°C) 

X24  *  Temperature,  950  mb  (10.1-15.5°C) 

x25  =  Temperature,  950  mb  (<10°C) 

X26  *  Temperature,  950  mb  (15.6-20°C) 
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A. 3  SAMPLE  FORECAST  PACKET  USED  AT  MOFFETT  NAS 


MOFFETT  STRATUS  FORMATION  FORECAST  AID  (SUMMER) 

The  following  computation  is  meant  as  a  forecast  aid  for  a  summer 
stratus  formation  forecast.  It  provides  guidance  for  the  stratus 
formation  forecast  which  should  then,  at  the  forecaster's  discretion, 
be  supplemented  by  other  information. 

The  computation  uses  the  0200GMT  observations  of  Moffett,  Alameda  am 
San  Francisco, and thaQOOO  GMT  Oakland  950mb  temperature  data.  Also  used  are 
"persistence"  ceilings  from  the  previous  day  at  Moffett.  Pive  forecast 
valid  times  (VT)  are  0600GMT,  0800GMT,  1000GMT,  1200GMT,  1400GMT,  which 
are  referred  to  as  computations  A  through  E. 

Erich  computation  produces  a  number  index,  which  falls  ideally  in  an 
interval  0.0  to  1.0,  however,  negative  numbers  and  numbers  larger  than 
1.0  do  occur.  The  forecast  index,  therefore,  extends  below  0.0  and 
above  1.0,  which  indicates  lower  and  higher  "probabilities"  respectively. 
Two  forecast  indices  are  computed  for  each  valid  time,  one  for  ceilings 
Z2000'  and  one  for  ceilings  41000'.* 

Following  the  computation  for  each  forecast  valid  time  is  a 
suggested  use  of  the  indices  and  the  interface  with  satellite  and 
synoptic  information.  Forecasters  may  develop  their  own  interface 
with  such  subjective  factors. 

The  percentages  appearing  on  the  right  of  the  "Forecast"  column 
are  the  verification  frequencies  (i.e.,  the  percentage  of  the  time  the 
forecast  was  incorrec^  in  the  developmental  data) . 

*  SYMBOLS 

Less  than 

Less  than  or  equal  to 
y -  Greater  than 
7  -  Greater  than  or  equal  to 
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The  following  data  will  be  needed. 

Moffett  (surface) 

Dew  point  depression  (0200GMT) 

Wind  Speed  (0200GMT) 

Previous  day's  ceiling  (0600,  0800,  1000,  1200,  6  1400GMT) 

San  Francisco  (SFO)  (surface) 

Wind  direction  (0200GMT) 

Dew  point  depression  (0200GMT) 

Alameda  (surface) 

Dew  point  depression 
Wind  speed 
Ceiling 

Oakland  (upper  air) 

Temperature  950mb  (OOOOGMT) 


(0200GMT) 

(0200GMT) 

(0200GMT) 
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SHEET  A 


MOFFETT  STRATUS  FORMATION  FORECAST  VALID  TIME:  0600  GMT 
Computation  for  ceiling  <  2000 ' 

If  the  following  parameters  fall  within  the  indicated  ranges,  add  or 
subtract  the  number  appearing  on  the  right. 

1.  Start  with  ,542 

2.  Moffett  dew.  pt.  dep,  (0.0-8.08F)  (0200  GMT)  .301  _ 

3.  Alameda  ceiling  (£2000?)  (0200GMT)  -.348  _____ _______ 

4.  Moffett  previous  days  ceiling:  0600GMT  (£2000')  -.155  _ 

Index  (TOTAL)  _ 

If  index  is:  Forecast 


£  .30 

Do  not  forecast  stratus  (<2000') 

2% 

.31 

to  .  40 

Do  not  forecast  stratus  unless  other 
indications  strongly  differ 

24% 

.41 

to  .80 

Need  additional  guidance  from  other 
indications 

43% 

>  .80 

Forecast  stratus  «2000') 

16% 

Enter  the  forecast  category  (jstratus/no  stratusl  on  the  final  forecast 
sheet  F. 
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Cent.  A 


MOFFETT  STRATUS  FORMATION  VALID  TIME:  0600GMT 
Computation  for:  Ceilings  ^1000 ' 

If  the  following  parameters  fall  within  the  indicated  ranges,  add  or 
subtract  the  numbers  appearing  on  the  right. 


1. 

Start  with 

2. 

Moffett  dew.  pt.  dep. 

(0 .0-8 . 0  °F)  (0200GMT) 

.183 

3. 

Alameda  ceiling  3£l000 

'  (0200GMT) 

-.253 

4. 

Moffett  previous  days 

ceiling:  0600GMT  (21000') 

-.141 

Index  (TOTAL) 


If  index  is: 

£.30 
.31  to  .40 

.41  to  .70 
*7.70 


Forecast 

Do  not  forecast  stratus  (<1000*)  4% 

Do  not  forecast  stratus  unless  other 
indicators  differ  28  & 

Need  additional  guidance  54% 

None  occurred 


Enter  forecast  on  sheet  F 


Sheet  B 


MOFFETT  STRATUS  FORMATION  FORECAST  VALID  TIME:  0800GMT 
Computation  for  ceilings  £2000' 

If  the  following  parameters  fall  within  the  indicated  ranges,  add  or 


subtract  the  number  appearing  on  the  right. 

1.  Start  with  .165 

2.  Moffett  dew.  pt.  dep.  (0.0-8.0°F)  (0200GMT)  .508  _ 

3.  Moffett  dew.  pt.  dep.  (8.1-10.0°F)  (0200GMT)  .220  _ 

4.  Alameda  ceilings  (X2000')  (0200GMT)  -.338  _ 

5.  Moffett  previous  days  ceilings  at  0800GMT  (22000')  -.129  _______ 

6.  Alameda  dew.  pt.  dep.  (0.0-8.0°F)  (0200GMT)  .110  _ 

7.  Alameda  wind  speed:  (12-15kts)  (0200GMT)  .131  _ 


Index  (TOTAL) 

If  index  is:  Forecast 


£.io 

Do  not  forecast  stratus  (<2000') 

8% 

.11 

to  .60 

Need  additional  guidance 

52% 

.61 

to  .70 

Forecast  stratus  (*.2000 ')  unless  other 
indications  strongly  differ 

24% 

y.70 

Forecast  stratus  U.2000') 

7% 

Enter  forecast  on  sheet  F. 


Cont.  B 


MOFFETT  STRATUS  FORMATION  FORECAST  VALID  TIME:  0800GMT 
Computation  for  ceiling  <1000' 

If  the  following  parameters  fall  within  the  indicated  ranges,  add  or 


subtract  the  numbers  appearing  on  the  right. 

1.  Start  with  _ .245 

2.  Moffett  dew.  pt.  dep.  (0.0-8.0°F)  (0200GMT)  .325  _ 

3.  Moffett  previous  days  ceiling  at  0800GMT  (Z1000 ' )-.  244  _ 

4.  SFO  wind  direction  (282°  -  8.0°F)  (0200GMT)  .107  _ 

5.  SFO  dew.  pt.  dep.  (0.0-8.0°F)  (0200GMT)  .096  _ 

6.  Oakland  950mb  temp.  (20-25 °F)  (0000GMT)  .075  _ 

7.  Alameda  ceiling  (Z10001)  (0200GMT)  -.146  _ 

8.  Moffett  wind  speed  (0-6kts)  (0200GMT)  -.104  _ 


If  index  is:  Forecast 


£..20 

Do  not  forecast  stratus  (£1000') 

3% 

.21 

to  .49 

Do  not  forecast  stratus  KIOOO') 
other  indicators  differs 

unless 

34% 

.50 

to  .60 

Need  addition  guidance 

49% 

.61 

to  .  80 

Forecast  stratus  (<£1000')  unless 
indicators  differ 

other 

35% 

7.80 

Forecast  stratus  (<.1000')  (O 

for  4) 

Enter  forecast  on  sheet  F. 
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Sheet  C 


MOFFETT  STRATUS  FORMATION  FORECAST  VALID  TIME:  1000GMT 
Computation  for  ceilings  <2000' 

If  the  following  parameters  fall  within  the  indicated  ranges,  add  or 
subtract  the  numbers  on  the  right. 


1. 

Start  with 

- 

2. 

SFO  dew.  pt.  dep, 

.  (0. 

0-8. 0°F)  (0200GMT) 

.212  _ 

3. 

Moffett  previous 

days 

ceiling  at  10 00 GMT  (22000 ') 

-.200  _ 

4. 

Moffett  dew.  pt. 

dep. 

(8 . 1-10 . 0°F)  (0200GMT) 

.185 

5. 

Moffett  dew.  pt. 

dep. 

(0 . 0-8 . 0°F)  (0200GMT) 

.694  _ 

Index  (TOTAL) 


If  index  is: 


Forecast 


£.30 


Do  not  forecast  stratus  (<2000')  11% 


.31  to  .40 


None  occurred 


.41  to  .50 


Do  not  forecast  unless  other  indicators 
differ  38% 


.51  to  .70 


Forecast  stratus  unless  other  indicators 
differ  39% 


7.70 


Forecast  stratus  (<£2000' ).  unless  other 
indicators  strongly  differ  26% 


Enter  forecast  on  sheet  F. 
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Cont.  C 


MOFFETT  STRATUS  FORMATION  FORECAST  VALID  TIME:  1000GMT 
Computation  for  ceilings  <1000* 

If  the  following  parameters  fall  within  the  indicated  ranges,  add  or 
subtract  the  numbers  on  the  right. 

1.  Start  with  - 

2.  Moffett  dew.  pt.  dep.  (0.0-8.0°F)  (0200GMT)  .352  - _ 

3.  Moffett  1000GMT  ceiling  previous  day  fcLOOO')  -.272  _ 

4.  SFO  dew.  pt.  dep.  <0.0-8.08F)  (0200GMT)  +.173  _ 


5.  Oakland  950rab  temp.  (10 . 1-15 . 5°F)  (0000GMT)  -.141 

6.  SFO  wind  direction  (237°-281°)  (0200GMT)  -.099 


Index 


(TOTAL) 


If  index  is:  Forecast 


4. 20 

Do  not  forecast  stratus  <1000' 

2% 

.21 

to  .  40 

Do  not  forecast  stratus  <1000'  unless 
other  indicators  strongly  differ 

20% 

.41 

to  .  50 

Do  not  forecast  stratus  <1000* 
other  indicators  differ 

unless 

24% 

.51 

to  .60 

Forecast  stratus  <1000*  unless 
indicators  differ 

other 

41% 

.61 

to  .  80 

Need  additional  guidance 

49% 

>.80 

Forecast  stratus  <1000'  unless 
indicators  differ 

other 

23% 

Enter  forecast  on  sheet  F. 
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SHEET  D 


MOFFETT  STRATUS  FORMATION  FORECAST  VALID  TIME:  1200GMT 
Computation  for  ceilings  <2000* 

If  the  following  parameters  fall  within  the  indicated  ranges,  add  or 
subtract  the  numbers  on  the  right. 

1.  Start  with  . 311 

2.  Moffett  dew.  pt.  dep.  (0.0-8.0®F)  (0200GMT)  .544  _ _ 

3.  Moffett  dew.  pt.  dep.  (>14.0°F)  (0200GMT)  -.121  _ 

4.  Moffett  dew.pt.  dep.  (8.1-10.0°F)  (0200GMT)  .248  _ 

5.  Moffett  previous  days  ceiling  at  1200GMT  (j>2000’>  -.130  _ 

6.  SFO  wind  direction  (282°-326°)  (0200GMT)  +.089  _ 

7.  SFO  dew.  pt.  dep.  (0.0-8.0°F)  (0200GMT)  .286  _ 

8.  Alameda  wind  speed  (12-15kts)  (0200GMT)  .116  _ 

Index  (TOTAL)  _ 

If  index  is:  Forecast 

<.30  Do  not  forecast  stratus  <2000*  8% 

.31  to  .40  Do  not  forecast  stratus  unless  other  16% 

indicators  strongly  differ 

.41  to  .49  Do  not  forecast  stratus  unless  other  26% 

indicators  differ 

.50  to  .90  Additional  guidance  needed  56% 

>.90  Forecast  stratus  (<  2000’)  15% 

Enter  forecast  on  sheet  F. 
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Cont.  D 


MOFFETT  STRATUS  FORMATION  FORECAST  VALID  TIME:  1200GMT 
Computation  for  ceiling  <1000* 

If  the  following  parameters  fall  within  the  indicated  ranges,  add  or 


subtract  the  number  on  the  right. 

1.  Start  with  _ •  325 

2.  Moffett  dew.  pt.  dep.  (0.0-8.0°F)  (0200GMT)  .414  _ 

3.  Moffett  previous  days  1200GMT  ceiling:  (£1000')  -.240  _ 

4.  SFO  dew.  pt.  dep.  (0.0-8.0°F)  (0200GMT)  .230  _ 

5.  SFO  wind  direction  (237°-281°)  (0200GMT)  -.095  _ 

6.  Alameda  dew.pt.  dep.  (>14.0°F) (0200GMT)  -.153  _ 

7.  Oakland  950mb  temp.  (15. 6-20. 0°F)  (0O00GMT)  -.121  _ 

8.  Oakland  950mb  temp.  «10.0°F)  (OOOOGMT)  -.248  _ 

9.  Oakland  950mb  temp.  (10.1-15.5)  (OOOOGMT)  -.257  _ 


Index  (TOTAL) 


If  index  is:  Forecast 

<.30  Do  not  forecast  stratus  <1000*  7% 

.31  to  .40  Do  not  forecast  stratus  <1000'  unless  other  36% 

indicators  differ 

.41  to  .70  Additional  guidance  needed  49% 

.71  to  .90  Forecast  stratus  <1000'  unless  other  33% 

indicators  differ 

>.90  Forecast  stratus  <1000'  15% 


Enter  forecast  on  sheet  F. 
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Sheet  E 


MOFFETT  STRATUS  FORECAST  VALID  TIME:  1400GMT 
Computation  for  ceilings  <2000' 

If  the  following  parameters  fall  within  the  indicated  ranges,  add  or 
subtract  the  number  on  the  right. 


1. 

Start  with 

.316 

2. 

SFO  dew.  pt.  dep. 

,  (0 . 0-8. 0°F)  ( 0200GMT) 

.186 

3. 

Moffett  previous 

days  ceilings:  1400GMT  (42000') 

-.162 

4. 

SFO  wind  direction  (237°-281°)  (0200GMT) 

-.124 

5. 

Moffett  dew.  pt. 

dep.  (8.0-10. 0°F)  (0200GMT) 

.132 

6. 

Moffett  dew.  pt. 

dep.  ( 0 . 0-8 . 0°F)  (0200GMT) 

.685 

7. 

Moffett  dew.  pt. 

dep.  (>14. 0°F)  (0200GMT) 

-.100 

Index  (TOTAL) 

If 

index  is: 

Forecast 

<  .20 

Do  not  forecast  stratus  <2000’ 

12% 

.21  to  .30 

Do  not  forecast  stratus  <2000*  unless 
other  indicators  strongly  differ 

26% 

.31  to  .60 

Additional  guidance  needed 

48% 

.61  to  .80 

Forecast  stratus  <2000'  unless 
indicators  strongly  differ 

other 

20% 

.81  to  .90 

Additional  guidance  needed 

52% 

>  .90 

Forecast  stratus  <2000' 

15% 

Enter  forecast  on  sheet  F. 
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Cont.  E 


MOFFETT  STRATUS  FORMATION  FORECAST  VALID  TIME:  1400GMT 


Computation  for  ceilings  <1000' 

If  the  following  parameters  fall  within  the  indicated  ranges,  add  or 


subtract  the  number  on  the  right. 

1.  Start  with  _ >476 

2.  Moffett  previous  days  ceiling,  1400GMT  t>1000')  -.257  _ 

3.  Moffett  dew.  pt.  dep.  (0.0-8.0°F)  (0200GMT)  .273  _ 

4.  SFO  wind  direction  (237°-281°)  (0200GMT)  -.071  _ 

5.  SFO  dew.  pt.  dep.  (0.0-8. 0°F)  (0200GMT)  .212  _ 

6.  Oakland  950mb  temp.  (10 . 1-15 . 5 °F)  (0000GMT)  -.167  _ _ 


Index  (TOTAL) 


If  index  is:  Forecast 


< 

.40 

Do  not  forecast  stratus  <1000' 

6% 

.41 

to 

.49 

Do  not  forecast  stratus  <1000*  unless  other 
indicators  differ 

33% 

.50 

to 

.89 

Additional  guidance  needed 

46% 

> 

.90 

Forecast  stratus  <1000' 

20% 

Enter  forecast  on  sheet  F. 


35 


Sheet  F 


VALID  TIMES  FORECAST  (<  2000')  FORECAST  (<1000*1 

0600GMT  _ 

0800GMT  _  _ 

1000GMT  _  _ 

1200GMT  _  _ 

1400GMT 


REMARKS : 

1.  Check  continuity  of  ceiling  forecast  throughout  the  period. 

2.  Check  for  ridiculous  results,  example:  Do  not  forecast  ceiling 
<1000',  without  also  forecasting  ceiling  <2000'  for  the  same 
valid  time. 
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APPENDIX  B 


LOGISTIC  REGRESSION  EQUATIONS  FOR  MOFFETT  NAS 


The  logistic  transf ormati on  is 


si 


6i  = 


1+e 


si 


1  ,  2 


where  9-j  is  the  probability  of  low  ceilings  ano  :  ^  =  1 1  - "  ■]  J  is 
the  probability  of  higher  ceilings. 

N 

s  .  =  I  A  .  .  X  . 

1  j-1  1J  J 

X  j  are  the  predictor  variables  and  the  A i j  are  the  constants 
found  by  using  regression  on  the  dependent  data  sample.  The  list 
of  candidate  predictors  and  forecast  equations  are  given  on 
succeeding  pages. 
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List  of  Variables 


=  East-west  component,  Moffett  wind  (kt) 

*2  =  North -south  component,  Moffett  wind  (kt) 

X3  =  Moffett  air  temperature  (°F) 

x4  =  Moffett  dew  point  depression  (°F) 

X5  =  Ceiling  (F.T.  minus  24  hrs)  Moffett  (Ft) 

x  6  =  Alameda  Ceiling  (Ft) 

X 7  =  East-west  component,  Alameda  wind  (kt) 

*8  =  North-south  component,  Alameda  wind  (kt) 

Xg  =  Alameda  air  temperature  (°F) 

x10  =  Alameda  dew  point  depression  (°F) 

X 1 1  =  San  Francisco  (SFO)  ceiling  (Ft) 

x12  =  East-west  component,  bFO  wind  (kt) 

X 1 3  =  North-south  component,  SFO  wind  (kt) 

x14  =  SFO  air  temperature  (°F) 

X 1 5  =  SFO  dew  point  depression  (°F) 

*lo  =  Marine  layer  thickness  (m) 

=  Marine  layer  strength  * 

x  18  =  700  mb  temperature  (°C) 

X 1 g  =  700  mb  wind  direction  (0-90°) 

x20  =  700  mb  wind  direction  (90-180°) 

Binary  variables 

X 2 1  =  700  mb  wind  direction  (  180-270°) 

x22  =  700  mb  wind  direction  (270-360°) 

X23  =  700  mb  wind  speed  (kt) 

x24  =  850  mb  temperature  (°C) 

^Defined  as  the  average  relative  humidity  from  the  surface  to  the 
marine  layer  top. 
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1 


Forecast  qua  t i on s 
Ceiling  j<  1 0 0 0  1 

Valid  Time 
2200  LT 

S  =  15. 319+. 134X]  -.096X2_. i67X3 

-.100X4-*634x10‘3X5--816X10‘3X6 
+  .068X7_.105X8-. 471X9+. 243X10 
+. 056 XI 2+ .416X14-. 615X15+. 010X18 
- .651X19+. 038X20+. 526X21-.  069X23 

2400  LT 

S  =  6.210+.113X1-.058X2_.296X4 
-.846X10-3X5-.  592X10-3X6+.  068X7 

-.170Xg  +  .  109XH-.  194X1  5+ -006X17 

0200  LT 

S  =  7.461+.083X r . 196X4_ . 5 10X10_3xb 
+  .086X7- .  343X9+.  19  IX 10+.  24  IX  14 

-*499X1 5-- 633X19+. 023X20+. 320X21 

+.013X24 

0400  LT 

S  =  6.093+.061Xr.196X4_#5ioxiO‘3x5 
+  .086X7- .343X9  +  .  191X10+.  24  1X14 
"• 499 xig-. 633X19  +  . 023X20+-  32 0X21 
+.013X24 

0600  LT 

S  =  6.137-.099X4_.7i7xi0"3x5 


+.079X7-. 071Xg -.056X12-. 045X13 
-. 1 99 Xi6_. 623X19-. 239X20+ -  656X21 

+ .007X24 


Forecast  Equations 
Cei 1 i ng  £2000 ' 

2200  LT 

S  =  .801-.  133X2-. 331X4-. 603X10"3X6 
1  12X8-.  356X15+  .  008X24 

2400  LT 

S  =  .  5  1 0+ .  1  50 X^  - .  0 7 3X2 - .  500X4 

-.  797X10-3X5+.  104X7-.  1  76X8 
-.24bxib_.  020X18-.  154X23+. 025X24 

0200  LT 

S  =  2.045+.126X1-.339X4_.645X10‘3X5 
- .151X8- .2 36X15 -.56 3X1 9- .124X20 
+  .522X21_. 077X23+. 012X24 

0400  LT 

5  =  .993+.091X] - • 283 X4 _ . 487 X i o " 3x5 
- -076X8+-8898X10_3Xli- .239X15 

- .  529Xig_ .088X20+. 605X21-  .069X23 
+ .014X24 

0600  LT 

S  =  -.100-.185X4_. 098X13-. 172X15 
+.012X24 
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APPENDIX  C 


LOGISTIC  REGRESSION  EQUATIONS  FOR  NORTH  ISLAND  NAS 

The  regression  equations  for  S  follow  from  the  list  of 
candidate  predictors  below.  Observation  time  is  1200  LT  unless 
otherwise  stated.  See  Appendix  B  for  definition  and  use  of  S. 

Li st  of  Predictors 

Xl  =  East-west  component,  North  Island  wind  (kt) 

x2  =  North-south  component,  North  Island  wind  (kt) 

X3  =  North  Island  air  temperature  (°F) 

x4  =  North  Island  dew  point  depression  (°F) 

X5  =  Ceiling  (FT  minus  24  hrs)  North  Island  (Ft) 

x6  =  East-west  component,  Miramar  wind  (kt) 

X j  =  North-south  component,  Miramar  wind  (kt) 

x8  =  Miramar  air  temperature  (°F) 

Xg  =  Miramar  dew  point  depression  (°F) 

x10  =  East-west  component,  San  Nicolas  wind  (kt) 

X^  =  North-south  component,  San  Nicolas  wind  (kt) 

x12  =  San  Nicolas  air  temperature  (°F) 

X13  =  San  Nicolas  dew  point  depression  (°F) 

x14  =  Marine  layer  height  (m) 

Xu  =  Marine  layer  strength  * 

x16  =  700  mb  temperature  (°C) 

Xu  =  700  mb  wind  direction  (degrees) 

♦Defined  as  the  average  relative  humidity  from  the  surface  to  the 
marine  layer  top. 
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Xig  =  700  mb  wind  speed  (kt) 

x19  =  850  mb  temperature  (°C) 

X20=  8^0  mb  wind  direction  (degrees) 

X21  =  850  mb  wind  speed  (kt) 

X22  =  Marine  layer  height  (Ft  minus  24  hr)  (m) 

x23  =  Surface  pressure  difference,  San  Nicolas 
North  Island  (mb) 

x24  =  Surface  pressure  difference,  Pt.Mugu  - 
North  Island  (mb) 

Forecast  Equations 
Ceil ing  <1000‘ 

Valid  Time 
1900  LT 

S  =  -4.713-.197X4+.513Xio-3X14 
+.022Xig 

0100  LT 

S  -  1. 624+. 080X^.091X3. .  107X4 

- . 005X15+. 026X19" -325X10_3X22 

0700  LT 

S  =  2 .783- . 106X4. .985X10"3x5 

-.00 3X15+. 015Xi 9+ .295X10-3x22 
+  . 046  X  24 
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Forecast  Equations 
Ceiling  £3500 1 

Valid  Time 
1900  LT 

S  =  10.58-.074X1-.122X3_.164X4 
-  .840X10-3x5+ . 002X15 -.043X18 

0100  LT 

S  =  11 .379- .141X3_ .08iX4 

-. 645X1 0-3X5+<013Xl6_>092Xl8 

0700  LT 

S  =  9.817  +  .O55X]-.14OX3_#732Xi0“^5 
+  .012X16-  .087X18+. 005X19 
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APPENDIX  D 


LOGISTIC  REGRESSION  EQUATIONS  FOR  MIRAMAR  NAS 

The  regression  equations  for  S  (see  Appendix  B  for  details) 
follow  from  the  list  of  predictors  below.  Observation  time  is 
1200  LT  unless  otherwise  stated. 

X1  =  East-west  component,  Miramar  wind  (kt) 

X2  =  North-south  component,  Miramar  wind  (kt) 

x3  =  Miramar  air  temperature  (°F) 

X4  =  Miramar  dew  po.nt  depression  (°F) 

x5  =  Ceiling  (FT  minus  24  hrs)  Miramar  (Ft) 

Xg  =  East-west  component,  North  Island  wind  (kt) 

x7  =  North-south  component.  North  Island  wind  (kt) 

X8  =  North  Island  air  temperature  (°F) 

x9  =  North  Island  dew  point  depression  (°F) 

Xio  =  East-west  component,  San  Nicolas  wind  (kt) 

XH  =  North-south  component,  San  Nicolas  wind  (kt) 

X}2  =  San  Nicolas  air  temperature  (°F) 

x13  =  San  Nicolas  dew  point  depression  (°F) 

X =  Marine  layer  height  (m) 

x15  =  Marine  layer  strength  * 

X i 5  =  700  mb  temperature  (°C) 

x17  =  700  mb  wind  direction  (degrees) 

♦Defined  as  the  average  relative  humidity  throughout  the 
marine  layer. 
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Xig  =  700  mb  wind  speed  (kt) 

*19  =  850  mb  temperature  (°C) 

X20  =  850  mb  wind  direction  (degrees) 

*21  =  850  mb  wind  speed  (kt) 

X22  =  Marine  layer  height  (Ft  minus  24  hr)  (m) 

*23  -  Surface  pressure  di fference ,  San  Nicolas 
North  Island  (mb) 

^24  =  Surface  pressure  difference.  Pt.Mugu  - 
North  Island  (mb) 


Val id  Time 


Forecast  Equations 
Ceiling  0000 ' 


1900  LT 

S  =  -1.608-.201X4>. 018X16 

+.020X19 

0100  LT 

S  =  I- 978-. 085X4_. 465X10-3X5 


+.043X18+. 016X19+. 053X24 

0700  LT 

S  =  - .553- . 07  7  X4_ .003X15 

- .003X15 pv+. 01 5X19+. 29 5X10- 3x22 

+.046X24 
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Forecast  Equations 
Ceiling  £3500 ‘ 

1900  LT 

S  =  8.532-.079Xr.073X3 
217X4-. 850xlC*3x5 

0100  LT 

S  =  7 .723- .07 IX y .07133- .074X4 

-.667X10-3X5+. 010X16-.  101X18-.  296X22 
+.043X24 

0700  LT 

S  =  10.273-.136X3_.745X10-3X5 
-.102X18+. 011X19 
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APPENDIX  F 


LOGISTIC  REGRESSION  EQUATIONS  FOR  SAN  NICOLAS  NAS 

The  regression  equations  for  S  (see  Appendix  B  for  details) 
follow  the  list  of  predictor  variables.  Observation  times  are 
1200  LT  unless  otherwise  stated. 

=  East-west  component,  San  Nicolas  wind  (kt) 

X2  =  North-south  component,  San  Nicolas  wind  (kt) 

x3  =  San  Nicolas  air  temperature  (°F) 

X4  =  San  Nicolas  dew  point  depression  (°F) 

*5  =  Ceiling  (FT  minus  24  hrs)  Miramar  (Ft) 

X6  =  East-west  component,  Miramar  wind  (kt) 

x7  =  North-south  component,  Miramar  wind  (kt) 

X8  =  Miramar  air  temperature  (°F) 

xg  =  Miramar  dew  point  depression  (°F) 

XiQ  =  East-west  component,  North  Island  wind  (kt) 

*11  =  North-south  component,  North  Island  wind  (kt) 

Xjg  =  North  Island  air  temperature  (°F) 

*13  =  North  Island  dew  point  depression  (°F) 

Xj4  =  Marine  layer  height  (m) 

*15  =  Ma"ine  layer  strength  * 

X}5  =  700  mb  temperature  (°C) 

*17  =  700  mb  wind  direction  (degrees) 

♦Defined  as  the  average  relative  humidity  throughout  the 
marine  layer. 
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X i g  =  700  mb  wind  speed  (kt) 

x19  =  850  mb  temperature  (°C) 

X20  =  250  mb  wind  direction  (degrees) 

*21  =  850  mb  wind  speed  (kt) 

X22  =  Marine  layer  height  (Ft  minus  24  hr)  (m) 

*23  =  Surface  pressure  difference,  San  Nicolas 
North  Island  (mb) 

x24  =  Surface  pressure  difference,  Pt.Mugu  - 
North  Island  (mb) 


Valid  Time 


Forecast  Equations 
Ceiling  _<1Q00 ' 


1900  LT 

S  =  8. 924-. 227X3+. 028X19 

0100  LT 

S  =  6. 282-. 137X3-. 001X5-. 001X14 
+.OI8X19 

0700  LT 

S  =  8 .434- .  1 6  3  X  3  _ . 00 1 X 1 4 
+.014X20+-432X23 
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Forecast  Equations 


Ceiling  £3500 1 

1900  LT 

S  =  14. 130-.152X16_. 249X8 

-.001X5+.  001X14+.  014X19 
+.047X24 

0100  LT 

S  =  13.242-.175X3_.001X5 

+  .009X16- .  167X18“ '588x33 

0700  LT 

S  =  8.970-  .093X3. .001X5 

-.118X18-  .048X24 
-. 217X4. ,850xl0“3x5 

0100  LT 

$  =  7.723-.071X  i-.07133. .074X4 

-.667X10-3X5+>01 ox 16-. 101X18-. 296X22 

+.043X24 

0700  LT 

S  =  10. 273-.  136X3-.  745X10-3X5 
- . 102X18+ .01 IX19 
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WASHINGTON,  DC  20235 


OFFICER  IN  CHARGE 
NAVOCEANCOMDET 
NAVAL  AIR  STATION 
MOFFETT  FIELD,  CA  94035 


COMMANDER  (2) 

NAVAIRSYSCOM 
ATTN:  LIBRARY  (AIR-00D4) 
WASHINGTON,  DC  20361 


USAFETAC/TS 
SCOTT  AFB,  IL  62225 


C  INCLANTT  l.T 

ATTN:  NSAP  SCIENCE  ADVISOR 

l '  'Df  N04E 

N.'CF  U..K  ,  VA  2351 1 


COMT  H !  P.DFLT 

ATTN:  NSAP  SCIENCE  ADVISOR 

CODE  N702/01T 

PEARL  HARBOR,  HI  96860 


COMSIXTHFLT/COMFAIRMED 
SCIENCE  ADVISORS  OFFICE  (032) 
FPO  NEW  YORK  09501 


OFFICER  IN  CHARGE 
OPTEVFOR,  SUNNYVALE 
NAVAL  AIR  STATION 
MOFFETT  FIELD,  CA  94035 


OFFICER  IN  CHARGE 
NAVOCEANCOMDET 
NAVAL  STATION 
FPO  SEATTLE  98791 


OFFICER  IN  CHARGE 
NAVOCEANCOMDET 
NAVAL  AIR  STATION 
SAN  DIEGO,  CA  92145 


COMMANDER 

NAVAIRSYSCOM  (AIR-330) 
WASH  I*:  3T0N,  DC  20361 


DIRECTOR  (12) 

DEFENSE  TECH.  INFORMATION 
CENTER,  CAMERON  STATION 
ALEXANDRIA,  VA  22314 


COMSECONOFLT 

ATTN:  NSAP  SCIENCE  ADVISOR 
FPO  NEW  YORK  09501 


COMSEVENTHFLT 

ATTN:  NSAP  SCIENCE  ADVISOR 
BOX  167 

FPO  SEATTLE  98762 


COMMANOER  NAVAL  AIR  FORCE 
U.S.  ATLANTIC  FLEET  (30F) 
ATTN:  NSAP  SCIENCE  ADVISOR 
NORFOLK,  VA  23511 


COMFLTAIR,  MEDITERRANEAN 
ATTN:  NSAP  SCIENCE  ADVISOR 
CODE  03A 

FPO  NEW  YORK  09521 


OFFICER  IN  CHARGE 
NAVOCEANCOMDET 
U.S.  NAVAL  AIR  FACILITY 
FPO  SEATTLE  98767 


OFFICER  IN  CHARGE 
NAVOCEANCOMDET 
NAVAL  AIR  STATION 
ALAMEDA,  CA  94501 


COMMANDER 
PACMISTESTCEN 
GEOPHYSICS  OFFICER 
PT.  MUGU,  CA  93042 


DIRECTOR 

TECHNIQUES  DEVELOPMENT  LAB 
GRAMAX  BLDG. 

8060  13TH  ST. 

SILVER  SPRING,  MD  20910 
WASHINGTON,  DC  20505 


52 


